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ABSTRACT 

We discuss the detection of large scale HI intensity fluctuations using a single dish 
approach with the ultimate objective of measuring the Baryonic Acoustic Oscillations 
(BAO) and constraining the properties of dark energy. To characterise the signal we 
present 3D power spectra, 2D angular power spectra for individual redshift slices, and 
also individual line-of-sight spectra, computed using the S 3 simulated HI catalogue 
which is based on the Millennium Simulation. We consider optimal instrument design 
and survey strategies for a single dish observation at low and high redshift for a fixed 
sensitivity. For a survey corresponding to an instrument with T sys = 50 K, 50 feed 
horns and 1 year of observations, we find that at low redshift (z w 0.3), a resolution 
of ~ 40 arcmin and a survey of ~ 5000 deg is close to optimal, whereas at higher 
redshift (z ss 0.9) a resolution of ~ 10 arcmin and ~ 500 deg 2 would be necessary - 
something which would be difficult to achieve cheaply using a single dish. Continuum 
foreground emission from the Galaxy and extragalactic radio sources are potentially 
a problem. In particular, we suggest that it could be that the dominant extragalactic 
foreground comes from the clustering of very weak sources. We assess its amplitude 
and discuss ways by which it might be mitigated. We then introduce our concept for a 
dedicated single dish telescope designed to detect BAO at low redshifts. It involves an 
under-illuminated static ~ 40 m dish and a ~ 60 element receiver array held ~ 90 m 
above the under-illuminated dish. Correlation receivers will be used with each main 
science beam referenced against an antenna pointing at one of the Celestial Poles for 
stability and control of systematics. We make sensitivity estimates for our proposed 
system and projections for the uncertainties on the power spectrum after 1 year of 
observations. We find that it is possible to measure the acoustic scale at z ss 0.3 with 
an accuracy ~ 2.4% and that w can be measured to an accuracy of 16%. 
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1 INTRODUCTION 

Precision cosmology is now a reality. Measurements of the 
Cosmic Microwave Background (CMB), Large-Scale Struc- 
ture (LSS) and Type la supernovae have led the way in defin- 
ing the standard cosmological model, one based on a flat 
Friedmann-Robertson- Walker Universe containing cold dark 
matter (CDM) and a cosmological constant, A, with densi- 
ties relative to the present day critical density f2 m « 0.27 
and S1a 



dark sector are usually geometrical (standard rulers or can- 
dles), or based on constraining the evolution of structure 
( Weinberg et al. [2012 1. 

The use of Baryonic Acoustic Oscillations (BAOs) to 
constrain w is now an established technique (Eiscnstci n et al.] 
1998] |Eisenstein|2003[ > . Originally applied to the Third Data 



Release (DR3) of the SPSS L uminous Red Galaxy (LRG) 
survey ( |Eisenstein et aL]|2005[ ), subsequent detections have 
been reported by a number of groups analysing data from 



0.73 (Riess et al. 



1998 



Perlmutter et al. 



1999 



2dF riCole et al. 2005b, SDSS DR7 (Percival et al. 20101, 



Eisenstein et al.||200l"] |Percival et al.||2002| |Komatsu eTat] WiggleZ (Bl ake et al.||2011[), 6dF ( |Beutler et al.||2011| ) and 



2011 1. The standard model has just 6 parameters, but these 
can been expanded to constrain, rule out, or detect more 
exotic possibilities. One particular case is an expanded dark 
energy sector, often parameterised by the equation of state 



parameter w — Pj p for the dark energy ( Copeland et al 
[2006]), or even the possibility of modifications to gravity 



( Clifton et al.pbll I. Observations which can constrain the 



BOSS ( |Anderson et ai.||2012| . All these have used redshift 
surveys performed in the optical waveband to select galaxies. 
The power spectrum of their three dimensional distribution 
is computed and, assuming that the power spectrum of the 
selected galaxies is only biased relative to the underlying 
dark matter distribution by some overall scale-independent 
constant, the acoustic scale can then be extracted. Present 
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sample sizes used for the detection of BAOs are typically 
~ 10 5 galaxies, but there are a number of projects planned 
to increase this to ~ 10 7 -10 8 using optical and near infra- 
red (NIR) observations which would lead to increased levels 
of statistical precision (Abbott et al. 2005 Hawken et al. 



2011||Laureijs et al.|2011| |Schlegel et al.|20li | 

It is possible to perform redshift surveys in the radio 
waveband, using 21cm radiation from neutral hydrogen (HI) 
to select galaxies which could provide improved confidence in 
the results from the optical/NIR. At present the largest sur- 
veys have only detected ~ 10 4 galaxies ( Meyer et al.|2004 |, 
limited by the low luminosity of the 21cm line emission 
( Furlanetto et al.|2006[ ) , but there are a number of proposed 
instruments with large collecting area capable of increasing 



this substantially in the next decade (Abdalla & Rawlings 
[2005} puffy et al.||2008"l |Abdalla et al.||2010[ ). By the time 
the Square Kilometre Array (SKA^j is online redshift sur- 
veys of galaxies selected in the radio should have comparable 
statistical strength to those made in the optical/NIR. The 
advantage of doing both radio and optical surveys is that 
they will be subject to completely different systematics. 

An intriguing possibility is the idea of using what has 
become known as 21cm intensity mapping. The fact that the 
21cm line is so weak means that it requires ~ 1 km 2 of col- 
lecting area to detect galaxies at z ~ 1 in an acceptable 
period of time - this was the original rationale for the SKA 
( Wilkinson || 199 1 1 . However, instruments with apertures of 
~ 100 m size have sufficient surface brightness sensitivity 
to detect HI at higher redshift, but they will preferentially 
detect objects of angular size comparable to their resolu- 



tion as pointed out in |Battye et al. (20041. These objects 
will be clusters with masses 10 -10 la M@ and other large 



scale structure. Peterson et al. ( 2006 1 and Loeb & Wyithe 



(2008) suggested that the full intensity field T(8, (j>, /) could 
be used to measure the power spectrum as a function of red- 
shift directly (as is done for the CMB intensity fluctuations) 
if the smooth continuum emission, e.g. that from the Galaxy, 
can be accurately subtracted and any systematic instrumen- 
tal artefacts, arising from the observational strategy, can be 
calibrated to sufficient accuracy. If these obvious difficulties 
can be mitigated, then this technique can be enormously 
powerful: the standard approach of detecting galaxies and 
then averaging them to probe large scale structure is inef- 
ficient in the sense that to detect the galaxies, at say 5a, 
one throws away all but a fraction of the detected 21cm 
emission. Although this is not a problem in the optical, in 
the radio where the signal is much weaker one cannot afford 
to take this hit. In order to measure the power spectrum 
one is not trying to detect the actual galaxies, although this 
can lead to a plethora of other legacy science, and all the 
detected emission is related to the LSS structure. |Chang| 
et al. ( 2008 I have shown that dark energy can be strongly 



constrained with a radio telescope significantly smaller, and 
cheaper, than the SKA. 

Interferometric techniques naturally deal with some of 
the issues related to foregrounds and instrumental system- 
atics. In particular, they can filter out more efficiently the 
relatively smooth and bright contaminating signals from 
the Galaxy and local signals such as the Sun, Moon and 
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ground. However, interferometers can be extremely expen- 
sive since they require electronics to perform the large num- 
ber of correlations. In this paper we will present our ideas 
on the design of a single dish telescope and an analysis of its 
performance (BINGO; BAO from Integrated Neutral Gas 
Observations). Attempts have already been made to detect 
the signal using presently operational, single-dish radio tele- 



scopes, ( Chang et el.||2010 | and recent surveys of the local 
universe (Pen et al 112008 1. It has been possible to make 



detections using cross correlation with already available op- 
tical redshift surveys, but not the detection of the crucial 
autocorrelation signal needed to probe cosmology using ra- 
dio observations alone. The reasons for this are unclear, but 
they are likely to be due to a combination of the obser- 
vations not been taken with this specific purpose in mind 
and the fact that the telescope and detector system are not 
sufficiently stable. Our philosophy in this paper will be to 
quantify the expected HI signal, quantify the contaminating 
foreground signals and then present our design for a system 
which will be sufficiently sensitive and stable for the BAO 
signal to be extracted. 



2 CHARACTERISATION OF THE SIGNAL 

In this section we will present some basic material on the 
21cm line emission in order to characterize various aspects 
of the signal expected in a HI intensity mapping experiment. 



2.1 Mean Temperature 

One can compute the emitted brightness temperature of the 
21cm line in a velocity width dv using 



Tcmitdf = 



hc A 3A 2 i 



(1) 



where k is Boltzmann's constant, T em i t is the HI brightness 
temperature, / cm it ~ 1420.4 MHz is the rest frame emission 
frequency for neutral hydrogen, nm is the number density 
of neutral hydrogen atoms, dl is the line of sight distance 
and A21 is the spontaneous emission coefficient of the 21cm 
transition. Equation (1) can be used to deduce the standard 
formula which relates the observed integrated flux from an 



object (typically a galaxy) with HI mass Mm ( Roberts 1974 \ 



M m = 



2.35 x 10 s M Q S obB dv 



l + z 



Jy kms~ 



Mpc 



(2) 



where d^(z) is the luminosity distance. 

This standard formula is for an isolated object. How- 
ever, on the very largest scales relevant to us here it is more 
convenient to work with the brightness temperature. The 
mean observed brightness temperature due to the average 
HI density in the Universe is 

„3 



T obs (2) 



^mit (z) _ f He 

l + z 



3^21 



( Phi (z) \ dl 



k 16/ c 2 mit Af H y \l + zjdv' 



(3) 

where Mh is the mass of the hydrogen atom and pm(z) is 
the average density of HI at redshift z. The first term is 
just a combination of fundamental constants and measured 
experimental parameters, the second is a function of redshift 
and the final, relatively unfamiliar, term relates the line of 
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sight distance to the recession velocity and is given by -ff^ 1 
in the local Universe, where Ho is the Hubble constant. 

In order to compute the last term in Equation j3b, let 
us consider the comoving volume element 



dV 



dzdQ H{z) 



"(*) = 



(1 + zf fJ , M2 fl + -j 

H(z) 



= c[d A (z)Y 



H{z) 



dz 





(4) 
(5) 



where d\(z) = r(z)/(l + z) is the angular diameter distance 
and . We can also write dV = d\dl dQ, and hence 



dl 



-(! + ; 



#(2) 

Therefore, using dv/c = dz/(l + z), we find that 



dv 



(X±*) 



(6) 



(7) 



Hence, (|3| becomes 



- / 3A 21 fc 3 \ (l + 0)Vx(g) 

Job S (,2J - TW72 — tjTZ\ ' W 



and substituting in values for the fundamental constants 
and re- writing pm{z) in terms of the density relative to 
the present day critical density, Qb.i(z) = 8irGpni{z)/(3Ho), 
gives 



2.45 x 10- 4 J E(z) ' 

where h = i^o/lOOkmsec" 1 Mpc -1 and E(z) = H(z)/H . 

Throughout we will assume that flmh = 2.45 x 1CP 4 
independent of redshift which is the value calculated from 



the local HI mass function measured by HIPASS ( Zwaan et 
20051). Both observ ational (|Prochaska fc Wolfe|2009| > and 



al. 



theoretical arguments ( Duffy et al.|20lT I appear to suggest 
that there is no more than a factor of 2 change in density 
of neutral hydrogen over the range of redshifts relevant to 
intensity mapping. In the rest of the paper we will drop the 
"obs" subscript and use T to denote the observed brightness 
temperature. 



2.2 3D power spectrum 

The 3D HI power spectrum can be written as 



[AT HI (fc,z)] 2 =f(zf[b(k,z)] 2 ^^ 2 



(k,z) 



(10) 



where b(k, z) is the bias and P c dm is the underlying dark 
matter power spectrum. We will assume, for the purposes 
of this discussion, that 6=1 independent of redshift, z, 
and scale k. 



Martin et al. 



(20121 find that this is a good 
assumption for scales > 10h~ Mpc implying that, within 
the constraints of their survey, that HI galaxies trace the 
large-scale structure of the Universe. In Fig. [I] we plot the 
predicted 3D linear power spectrum of 21 cm fluctuations 
for z — 0.28 which corresponds to the case (A) discussed 
in section [3] - the spectrum at other values of z is virtu- 
ally indistinguishable on this logarithmic scale except for 
the overall growth factor. At low k, large angular scales, the 
signal has a spectrum Thi oc fc 2 , which turns over around 



Figure 1. The 3D HI power spectrum at z = 0.28. BAOs are 
just about visible in the range 0.01 < fc/Mpc -1 < 0.2, but they 
are less than 10% of the total signal. For k < 0.01 Mpc -1 the 
temperature spectrum is oc k 2 corresponding to P(k) oc fc, that 
is, scale invariant. 



fc ~ 0.1 Mpc -1 . BAOs are less than 10% of the signal in 
the range 0.01 < fc/Mpc" 1 < 0.2 and the power spectrum 
is significantly affected by the non-linear regime for fc > k n \ 
which is « 0.2 Mpc -1 at z — 0. The wavenumber fc can be 
converted to angular multipole using £ ~ fcr(z) where r(z) 
is the coordinate distance and this can in turn be converted 
to an approximate angle scale using 9 m 180 deg/£. The 
non-linear scale corresponds to £ B \ w 220 and 9 n i w 0.8 deg 
at low redshift [z ~ 0.3). Observations which probe signifi- 
cantly higher resolution than this will require the modelling 
of non-linear effects on the power spectrum. Note that this 
is somewhat different at higher redshifts. 



2.3 2D Angular power spectrum 

One can construct a 2D angular power spectrum for the 
HI intensity over some frequency range which can then be 
treated in using techniques designed to study the CMB 
anisotropies. We start from the 3D quantity ST(r(z)n, z) 
which can be computed directly from (JsJ) by replacing pm 
with Spiii- This makes the assumption that d£/dv and 1 + z 
take the values in an unperturbed universe, that is, we ignore 
the effects of peculiar velocities, which will lead to redshift 
space distortions that could be important on small scales, 
and the Sachs- Wolfe effects which could be important on 



large scales ( |Challinor fc Lewis 2011 Bonvin & Durrer 
20111. Its projection on the sky ST(h) is defined using a 
normalized projection kernel W(z), that is 



<ST(n) = / dzW{z)ST(r(z)h,z) 

dzW(z)f{z)5 H i(r{z)n,z). 
We will use a top-hat projection kernal 

- - if Z m in ^ Z ^ Z n 



(11) 



W(z) 



otherwise, 



but the concept is more general and any weighting function 
could be used. The maximum and minimum redshifts of the 
bin are related to the minimum and maximum frequencies 
of the bin. 
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Fourier transforming and then expanding the Fourier 
modes in spherical harmonics, we have 



5T(r(z)n,z) 



T(z) 



d s k - 



(2tt)= 



5 H i(k, z)e 



ir (z)nk 



d k - 

j^S H i(k,z)j e (kr(z)) 



l,m 

x F/ m (k)Vi m (n) , (12) 
where je is the spherical Bessel function. We then have 

ST(h) = 4tt^Y [ dzW(z)f(z) 



I ,m 

d 3 k 
(2^r)3 



fe,(M)i*(M*))*ta,(k)*im(A) , (13) 



which gives 



a lm = 4vrr / dzW(z)T{z) 



d*k f 



(27T) 



-5jj/(k, «)j^(fcr(2))y/ m (k) 
(14) 

The angular power spectrum Ci is found by taking the 
ensemble average 



(15) 



We need to use the orthonormality condition for the spher- 
ical harmonics, statistical isotropy and the definition of the 
power spectrum 

(5 HI {k,z)S H i(k',z')} = (27vfS(k--k')b 2 P cdm {k)D(z)D(z'). 

(16) 

Here, D(z) is the growth factor for dark matter perturba- 
tions defined such that D(0) = 1. Then we find that 



2 -k- J dzW(z)T(z)D(z) J dz'W(z')T(z')D(z') 

x f k 2 dkP cdln (k)j e (kr(z))j e (kr(z')). (17) 



For large £, we can use the Limber approximation ( Lim- 
ber||1953 Afshordi fc Loverde||2008 \ to perform the k inte- 
gral, that is, 



k 2 dkP cdm (k)ji(kr)ji {kr') 



— Pcd 



l + l/2\ S(r-r') 



(18) 
(19) 



Writing cdz — HoE(z')dr' , we can use the 5-function to 
perform the dr' integral and we deduce 



c t 



dzE(z) 



W(z)T(z)D(z) 



r(z) 



(20) 

We can proceed to calculate the 21cm angular power 



spectrum using (20 1. The matter power spectrum P c dm(k) 
today (z = 0) can be computed from 



Pcdm(fe) = Ak nS T 2 (k). 



(21) 



For the transfer function T(k) we use the fitting formu- 
lae given by Eisenstein & Hu (19981, assuming Q m = 0.27, 
SIa = 1 — fi m , h = 0.71 and baryon fraction /baryon = 0.17. 
We take the primordial spectral index to be ns = 0.96 and 



use the approximate formula by Bunn & White (19971 for 



|925-975| MHz 




Figure 2. The HI angular power spectrum for the same central 
frequency, 950 MHz, but for two different ranges 925-975 MHz 
and 800-1100 MHz. The signal amplitude increases with decreas- 
ing the frequency range and the BAOs become more prominent. 
As the frequency range increases the size of the region probed 
becomes larger implying that the fluctuations will decrease. 



the power spectrum amplitude. For this choice of param- 
eters, the sound horizon s ~ 152 Mpc and the BAO scale 
k A = 2tt/s w 0.041 Mpc" 1 . The growth function D(z) is 
calculated using the fitting formula from |Carroll, Press fc] 



Turner 



In 



(1991). 



we plot the HI angular power spectra for the 
same central frequency 950 MHz (corresponding to a cen- 
tral redshift z c ~ 0.5) but for two different bin widths 
A/ = 300 MHz and A/ = 50 MHz. We see that as a re- 
sult of decreasing the bin size, the signal amplitude in- 
creases and the BAO wiggles become more prominent. It 
is clear that there exists an optimum frequency bin size 
for which a detection will be possible: the signal-to-noise 
level would increase oc y/Af if the signal level was inde- 
pendent of frequency, but the level of fluctuations in the 
HI is reduced as the box enclosed by the angular beam-size 
and frequency range are increased. In Fig. [3] we plot the ra- 
tio Ci/CV, smoo th, where Ci, Bmoo th is the no-wiggles angular 



power spectrum ( Eisenstein & Hu 1998 1 for the same central 



frequency 950 MHz but for the two frequency ranges used in 
Fig. [3] In addition we have plotted the ratio -P(fc)/P smoot h(k)> 
where P(k) — P c dm(k). Note that the ratio C^/CV iSmoot h 
cannot be bigger than P(fc)/P smoot h(k)- However, the two 
are in close agreement for the Af = 50 MHz bin which indi- 
cates that this is the optimal frequency bin size for tracing 
the BAOs. 



2.4 Estimated signal from S 3 simulations 



Obreschkow et al. ( 2009a b ) present semi-analytic simu- 
lations of HI emission based on catalogues of galaxies 
whose properties are evolved from the Millennium or Milli- 



Millennium dark matter simulations (Springel et al. 20051 



and we will use this to compute simulated spectra as would 
be observed by a single dish. We adopt the smaller of 
the two resultant mock sky fields (corresponding to the 
Milli-Millennium simulation, comoving diameter Sb ox = 
62.5 h~ Mpc) for the characterisation of the expected HI 
signal which, as we will discuss later, restricts the highest 
redshift to which the spectrum will be complete. We extract 
galaxy properties from the simulated catalogue - position co- 
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1800-1 100| MHz 




k |Mpc" | 



Figure 3. In the top plot we present the ratio Ci/Ci smooth f° r 
the same central frequency (950 MHz) but for different frequency 
ranges. In the bottom we present the ratio P(k) / P(k) SYnoot ^. The 
925 — 975 MHz curve in the upper plot is very similar to that in 
the lower one implying that a frequency range of 50 MHz will be 
optimal to detect the BAO signal at this redshift. 



ordinates, redshift, HI mass, integrated HI flux and Hi-line 
width - from which we calculate the HI brightness tempera- 
ture contributed by galaxies in the frequency range A/, as 
it would be observed by a telescope beam which is a top 
hat of diameter #fwhm independent of frequency. Using a 
diffraction limited beam for which #fwhm oc 1 + z would 
result in a higher signal at lower frequencies. 

In Fig. [4] we present the expected HI temperature fluc- 
tuations between frequencies of 800 and 1200 MHz, which is 
part of the frequency range we will be interested in, for a 
range of beam sizes, Afwhm and frequency bandwidths, A/. 
These plots provide considerable insight into the nature of 
the signal that we will be trying to detect. It appears that the 
average signal from equation Q, represented by the dashed 
line, is compatible with the spectra from the Millennium 
Simulation. Larger values of 6fwhm and A/ mean that the 
volume of the universe enclosed within a single cell is larger 
and hence the spectra have only very small deviations from 
the average, whereas for smaller values there are significant 
deviations as large as 400 yX. illustrating that the signal on 
these scales can be much higher than the average. 

There are limitations to this approach: our spectra are 
not accurate out to arbitrarily high redshift where the cone 
diameter exceeds Sbox and hence simulated catalogue is miss- 
ing galaxies, that is, it is not complete. The maximum red- 
shift for which the simulated catalogue is complete is de- 
termined by the beam size, with smaller beam sizes being 





(a) A/ = 1 MHz, $fwhm = 60arcmin 




Frequency |MHz] 

(b) A/ = IOMHzJpwhm = 60arcmin 





wmm 



1000 

Frequency [MHz] 



(c) A/ = 1 MHz, e FWHM = 20arcmin 




(d) A/ = 10 MHz, PWH M = 20 arcmin 

Figure 4. The simulated HI signal for different telescope beam 
resolution and frequency bin size combinations, Sfwhm an d A/ 
respectively. A smaller frequency bin or telescope beam size re- 
sults in larger fluctuations, as expected, since going to smaller 
values of SpwHM and A/ means that one is probing a smaller 
volume of the Universe and fluctuations in the HI density are ex- 
pected to be larger on smaller scales, as they are in the cold dark 
matter. In each case the dashed line is the average temperature 
computed using |9jl. The objective of 21cm intensity mapping ex- 
periments is to measure the difference between the signal and this 
average. 
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Figure 5. The number of galaxies contributing to the HI signal 
in each bin as a function of redshift. We have used Sfwhm = 
20 arcmin and A/ = 1 MHz for this particular example. Note 
that the signal presented in Fig. [4] varies only very weakly with 
redshift yet the number of galaxies contributing to each bin in- 
creases dramatically. At low redshifts, we see that the signal is 
dominated by individual galaxies, at intermediate redshifts it is 
dominated by galaxy clusters and groups as suggested in | Battye] 
|et al.| | |2004| and at high redshifts the signal comprises an unre- 
solved background signal, that is, there arc a very large number 
of galaxies contributing to each bin in redshift. 



accurate to higher redshifts. For a beam size of 1°, the spec- 
trum will be accurate up to around 600 MHz. 

At the other end of the frequency range, which is not 
shown in Figj4] there are large fluctuations around the 
mean that can be attributed to the effects of shot noise. In 
Fig. [5] we plot the number of galaxies contributing to each 
bin of the spectrum defined by #fwhm = 20 arcmin and 
A/ = 1MHz, that is, Fig. |4jc), as a function of redshift. 
At low redshifts (high frequencies) the number of galaxies 
contributing to each bin is small, and the signal cannot be 
thought of as an unresolved background. If N is the number 
of contributing galaxies, the deviation due to shot noise is 
a 2 oc l/Af and consequently, we see larger fluctuations. This 
is exactly the regime probed by surveys such as HIPASS 
and ALFALFA which detect individual galaxies for z < 0.05 
with slightly better - but not substantially better - resolution 
than used in this plot. As one moves to higher redshifts a low 
resolution telescope will preferentially detect objects which 
have an angular size which is similar to that of the beam. 
At intermediate redshifts this will be the clusters discussed 



in Battye et al. (20041 and these can be clearly seen as the 
well-defined spikes in the spectrum around z = 0.3 — 0.4. 
At higher redshifts even clusters are beam-diluted; beyond 
z — 0.5 there are always a few galaxies contributing to each 
bin - essentially there is an unresolved background. 

We can quantify the HI temperature fluctuations as a 
function of #fwhm and A/ by calculating their root mean 
squared (r.m.s.) deviation about the mean over the fre- 
quency range. As shown in Fig. [6] the r.m.s. decreases with 
increasing Afwhm and A/. This is exactly what is expected 
since as #fwhm and A/ increase the region of the universe 
enclosed is larger and the observed temperature will become 
closer to the average. We note that typically the fluctuations 
are significantly bigger than one for a substantial part of the 
parameter space. 




Figure 6. The r.m.s. temperature fluctuations in the HI signal 
as a function of $fwhm (solid line) and A/ (dashed line). When 
we vary #fwhm we nx A/ = 1 MHz and in the case where A/ is 
varied we fix #fwhm = 20 arcmin. 



3 OPTIMIZING THE PERFORMANCE AND 
SCIENCE POTENTIAL OF A SINGLE DISH 
EXPERIMENT 

In this section we analyse the expected performance of a 
single dish telescope in two frequency ranges either side of 
the mobile phone band around 900-960 MHz: (A) / = 960- 
1260 MHz (0.13 < z < 0.48) and (B) / = 600 - 900 MHz 
(0.58 < z < 1.37). We will compute projected errorbars on 
the resulting power spectra based on the methods presented 
in Seo et al. (2010 \ using the 3D power spectrum ([!]). In each 
case. (A) and (B), we will optimize the area covered, Q surt 



and the resolution of the telecope, #fwhm oc l/D, where D 
is the illuminated area of the dish, in order to give the best 
possible measurement of the acoustic scale. 

The standard expression for the projected error on a 
power spectrum measurement averaged over a radial bin in 
fc-space of width AA; is ( Feldman, Kaiser fc Peacock ||199~ 
Seo et al.||2010 ) 



OP = / (277)3 1 

P V V SUI 4nk 2 Ak 



J2 17 
O nix ' pix 



1+ _ 

[T(z)] 2 W(k) 2 P 



(22) 



where we have neglected the contribution from shot noise. In 
order for the analysis to be simplified, it is necessary for us to 
make the errorbars on these so-called bandpower measure- 
ments uncorrelated. This requires that Ak > Tr/(r(z)9 nl i n ) 
where we assume that the observed sky area is a square 
patch with f2 sur = 6> m i n x 6> max and # max > m i n . In order to 
be able to detect the acoustic scale we need Ak to be less 
than a few times k\ and this can impose important limi- 
tations on the design of drift scan surveys where typically 
#max 3> #mm- V SU r is the survey volume given by 

V sm = Q sm f dz-^, (23) 



dzdQ. 



where 



dV 



dzdQ, H E{z) 



(24) 



The window function W(k) models the instrument's an- 
gular and frequency response. Because our frequency resolu- 
tion is relatively high, we can ignore the instrument response 
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Figure 7. Optimization of the measurement of &a for case (A) - 
observations over the frequency range 960 — 1260 MHz. We have 
computed the projected uncertainty on the measurement of &a as 
a function of Af2 sur and the beam size SpwfjM- In the top panel 
the contours are lines of constant <5&a with the numbers on the 
contours signifying the percentage fractional error: IOO^a/^A- I n 
the middle panel we present a ID slice as function of $fwhm 
with Af7 sur = 2000 deg 2 and in the bottom panel a ID slice as a 
function of f2 SU r with fpwHM = 40 arcmin. 




0FWHM [deg] 




a m [deg J ] 

Figure 8. Optimization of the measurement of k\ for case (B) 
- observations over the frequency range 600 — 900 MHz. We have 
computed the projected error on the measurement of k\ as a 
function of Af! sur and the beam size %whm- I n the top panel the 
contours are lines of constant Skj^. The numbers on the contours 
are defined in the same way as in Fig. [7] In the middle panel we 
present a ID slice as a function of Sfwhm with Af! sur = 500 deg 2 
and in the bottom panel a ID slice as a function of Q B ur with 
Opwhm = 10 arcmin. 



function in the radial direction. However, we have a finite 
angular resolution, which imposes a window function 

2" 



W(k) = exp 



1 



k r(z) 



'FWHM 



(25) 



The beam area is given by f2 p i x oc #p WHM and the pixel 
volume is given by 

z+Az/2 



dz 



z-Az/2 



dV 
dzdQ. 



(26) 



where Az is the redshift range corresponding to A/, the 
nominal frequency channel width, which we will take to be 
1MHz. 

The pixel noise for a perfect, dual-polarisation, switch 
beam receiver measuring is given by 



(27) 



where T sys is the overall system temperature including con- 
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Figure 9. Projected errors on the power spectrum (divided by 
a smooth power spectrum) expected for the survey described by 
the parameters used in the text and Table [I] which is close the 
optimal value for Fig[7] that is, case (A) from Table [l] We have 
used AA; = 0.016 Mpc -1 . The projected errors would lead to a 
measurement of the acoustic scale with a percentage fractional 
error of 2.4%. Note that the points do not lie on the line since 
they are band powers and hence are integrals over neighbouring 
values of k. 




Figure 10. Projected errors on the power spectrum (divided by 
a smooth power spectrum) expected for the survey described by 
the parameters used in the text and table [T] which is close the 
optimal value for Fig|8] that is case (B) from table [l] We have 
used Afc = 0.006 Mpc -1 . The projected errors would lead to a 
measurement of the acoustic scale with a percentage fractional 
error of 1.5%. 



tributions from the sky, ground and receiver, which we will 
assume to be 50 K - a value that is conservative for a well 
designed system based on room temperature low noise am- 
plifiers (LNAs) - and t plx is the time spent observing one 
pixel. This is related to the total observing time, t b a , by 



tobs — 



(28) 



where uf is the number of feedhorns at the focus of the 
telescope. In what follows we will keep t b s and 71f fixed 
to be 1 year, that is a full year of on-source integration, 
something which may take around 2 years of elapse time to 
achieve, and 50, respectively. In order to estimate the error 
on the measurement of k A we will follow Blake & Glazcbrook 
( 2003 1 and fit a decaying sinusoidal function to the projected 



data 



P(k) 



= 1 + Ak exp 



k 



0.1 ft Mpc" 



(2nk\ 



V fc A / 



where the two parameters are the acoustic scale, k\, and the 
overall amplitude, A, which we will ultimately marginalize 
over. 

The parameter A; a is a function of the cosmological pa- 
rameters and in particular the matter density measured rel- 
ative to the critical density, il m , and the equation of state 
parameter of the dark energy, w — Pdc/Pdo = Wo + w\z. 
The angle-averaged power spectrum allows one to measure 
the isotropic BAO scale k\ = (fei^y) 1 '"' 3 m terms of the 
wavenumbers parallel to, kit, and perpendicular to, k±, the 
line-of-sight. An efficient and quick method to do this is to 
use the approximations ( Seo fc Eisenstein|20 03 Glazebrook 
fc Blake|200"5 l 



fei 



H( 



j r 

-fill 



. kt 



<7,\ 



Z) , ref 

k i 



(30) 



in terms of reference values denoted by "ref used for nor- 
malization purposes. Our reference cosmology is our input 
cosmology, with w — —1 = const, fi m = 0.27 and h = 0.71. 
This allows us to define the angle-averaged distance to red- 
shift 2 is given by ( |Beutler et al.|2011| > 



d v (z) = d A (z) 



. cz(l + z) 
H(z) 



1/3 



(31) 



which takes into account the angular diameter distance, 
<1a{z), and the radial distance c/ 'HqE(z). 

In Figs. [7] and [8] we present projected percentage errors 
on the acoustic scale, 100<5fcA/fcA as a function of the sur- 
vey area, O sur , and the beam size, #fwhm for the redshift 
ranges given by cases A and B. First it is clear that in each 
case there are a wide range of survey parameters for which 
accurate measurements of A?a can be made. It will always 
be the case that reducing 6>fwhm will lead to a better mea- 
surement of k\. However, beyond some point one gets into 
a situation where there is a very significant extra cost for 
very little extra return. 

In case A for Q B ur > 2000 deg 2 the contours of constant 
Skx/kA are almost vertical. This implies that for a fixed 
value of #fwhm , so long as one can achieve f2 sur > 2000 deg 2 , 
optimisation of the survey area can only yield minor im- 
provements. This is illustrated in the bottom panel of Fig. [7] 
where we fixed #fwhm- The optimum f2 sur is between 5000 
and 6000deg 2 where Skx/kA ~ 0.0225, but for example if 
S7 SUI . = 2000 deg 2 , as we will suggest for the BINGO instru- 
ment discussed in Section [5] SkA/kx ~ 0.024, a difference 
which is within the uncertainties of the present calculations. 
The exact position of the optimum is governed by a balance 
between cosmic variance and thermal noise terms in equa- 
tion (22 1; typically, it is where the two are approximately 



(29) 



equal. In the middle panel of Fig. [7] we fix Q BUT = 2000 deg 2 
and vary only #fwhm- There is a significant improvement 
in 5kA/k\ as one goes from #fwhm = 2 deg to around 
40 arcmin but beyond that point there is little improve- 
ment. The reason for this is made apparent in Fig. [9] where 
we present the projected errors on P(fc)/P smoot h(A;), where 
^smooth(fe) is a smooth power spectrum chosen to highlight 
the BAOs. We see that as one goes to higher values of k, cor- 
responding to higher angular resolution and lower values of 
Sfwhm, the overall amplitude of the BAO signal decreases 
and therefore it is clear that the discriminatory power, in 
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Figure 11. Projected constraints on the Hubble diagram for the 
volume averaged distance, dy(z), for cases A and B described by 
the parameters in table [T] Included also are the actual measure- 
ments made by 6dF, SDSS-II, BOSS and WiggleZ. The top panel 
is the absolute value and the bottom panel are the residuals from 
the fiducial model. On the bottom panel the shaded region indi- 
cates the range of dy allowed by the lc constraint C m /i 2 from 
WMAP7 | |Komatsu et al.|2011^ . The dotted line is the prediction 
for wo = —0.84 (which could in principle be ruled out by case A) 
and the dashed line is too = —0.93 (by case B). 



terms of measuring BAO^j comes from the first three oscil- 
lations in the power spectrum. Although measuring beyond 
this point will improve the constraint, the improvement does 
not justify the substantial extra cost involved; going from 
Sfwhm = 40 arcmin to 20 arcmin would require the telescope 
diameter to increase by a factor of 2 and the cost by a factor 
of somewhere between 5 and 10. For the subsequent discus- 
sion of case A, we will concentrate on Afwhm = 40 arcmin 
and fi sur = 2000 deg 2 and the relevant survey parameters 
are presented in Table [l] . This is not exactly the optimum 
point in the #fwhm — f2 aur plane, but it is sufficiently close 
to the optimum. #fwhm = 40 arcmin can be achieved with 
an illuminated diameter of around 25 m at the relevant fre- 
quencies and r2 sur = 2000 deg 2 can be achieved with a drift 
scan strip of around 10 deg at moderate latitudes which is 
compatible with our choice of Ak. We note that the precise 
values of the optimum depend on the overall survey speed 
which is governed by t b B nF /T s 2 ya . Typically an increase in 
the overall survey speed will move the optimum survey area 
to a larger value. 

The situation for case (B) is qualitatively similar, but 
the exact values are somewhat different. Being centred at a 



much higher redshift, z ~ 0.9, means that the angular size of 
the BAO scale is much smaller and hence higher resolution is 
necessary to resolve the BAO scale and a smaller area needs 
to be covered to achieve the same balance between thermal 
noise and cosmic variance. We find that something close to 
optimal can be achieved for 6>fwhm = 10 arcmin and Q SUI = 
500 deg 2 . The projected errors for such a configuration are 
displayed in Fig. 10 which would lead to SkA/kA ~ 0.015. In 
order to achieve this, one would need a single dish telescope 
with an illuminated diameter of ~ 100 m which would have 
to be steerable, at least in some way. It would be relatively 
expensive to build a dedicated facility to perform such a 
survey, although, of course, such facilities do already exist, 
albeit with very significant claims on their time from other 
science programs. Perhaps an interferometric solution could 
be used to achieve this optimum: we find that an array of 
~ 250 antennas with a diameter of 4 m spread over a re- 
gion of diameter ~ 100 m, having the same value T sya , could 
achieve the same level of sensitivity. More generally an inter- 
ferometer comprising rid antennas of diameter D spread over 
a region of diameter d would be equivalent to a single dish 
of diameter d with a focal plane array with np feedhorns, if 
rid = d v / n^/D so long as np < nd, that is the filling factor 
of the interferometer is < 1. 

In order to get an impression of the science reach of 
the proposed observations we have done two things. First, 
in Fig. [Tl]we have plotted our projected constraints on the 
angle-averaged distance on the BAO Hubble diagram - both 
absolute and residuals - along with the state-of-the-art mea- 
surements from 6dFGS, SDSS-II, BOSS and WiggleZ. It is 
clear that the measurements that would be possible from 
instruments capable of achieving the design parameters of 
cases A and B would be competitive with the present crop 
of observations and those likely to be done in the next few 
years. They would also be complementary in terms of red- 
shift range and add to the possibility of measuring the cos- 
mic distance scale solely from BAOs. We note that it should 
be possible, in fact, to make measurements of the acoustic 
scale over the full range of redshifts probed by the intensity 
mapping surveys and it would be possible to place a num- 
ber of data points on the Hubble diagram, albeit with larger 
errorbars. In order to compare with the present surveys we 
decided to present a single data point. 

In addition, we have made simple estimates of the 
projected errors of the dark energy parameters assuming 
that all the other cosmological parameters, including Q m 
are fixed. This will be a good approximation in the era 
of high precision CMB measurements made by, for exam- 
ple, the Planck Satellite]^] . For constant w, the measure- 
ment of SkA/kA ~ 0.024 possible in case A would lead to 
Sw/w — 0.16, whereas for case B, where SkA/kA = 0.015 
would lead to Sw/w = 0.07. For models with non-constant 
w, we have performed a simple likelihood analysis for the 
two-dimensional parameter space wo — w i and the results 
are presented in Fig. [12] for cases A and B and a combina- 
tion of the two. We see that since we have assumed a single 
redshift measurement of the acoustic scale, there is an exact 



2 We note that this might be different if we had chosen another 
science goal such as the measurement of neutrino masses. 



3 We note that varying O m h 2 modifies dy(z) and the present 
limits on this parameter combination are included in Fig. |11| 
These will be substantially improved in the next couple of years. 
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Table 1. Table of the survey parameters used in Figs, |9]and |10| In both cases we have chosen the values of 0FWHM and Q SU r so as to 
be close the optimal values as described in the text. Included also are the projected measurements of the acoustic scale and estimates of 
the error w (assumed constant) when all other cosmological parameters are known. 



Property 


Low redshift (A) 


High redshift (B) 


^centre 


0.28 


0.9 


^(^centre ) 


1100 Mpc 


3000 Mpc 


^sur 


2000 deg 2 


500 deg 2 


Vsur 


1.2 Gpc 3 


3.1 Gpc 3 


#FWHM 


40 arcmin 


10 arcmin 




810 Mpc 3 


590 Mpc 3 


r.m.s. noise level (A/ = 1 MHz) 


85 fiK 


170 fiK 




0.024 


0.015 


Sw/vj 


0.16 


0.07 




Wr, 

Figure 12. Projected uncertainties on the dark energy pa- 
rameters wo and u>i for case A with $fwhm = 2/3 deg and 
O sur = 2000 deg 2 (top), case B with Sfwhm = 10 arcmin and 
f2sur = 500 deg 2 (middle) and a combination of the two (bottom) 
for the case where all other cosmological parameters are known. 
We have assumed that cases A and B lead to a measurement 
of the acoustic scale at the central redshift of the survey which 
means that there is an absolute parameter degeneracy between 
the two parameters, but the degenerate line is different at the two 
redshifts. Combination of the two removes the degeneracy. 



degeneracy between the two parameters, but the direction 
of the degeneracy is very different at the two different red- 
shifts. When they are combined the degeneracy is broken 
and it is possible to estimate the two parameters separately. 



4 FOREGROUND CONTAMINATION 
4.1 Overview 

To detect BAO in the faint HI signal, we will need to remove 
the contributions from much brighter foreground emission. 
The observed brightness temperature of the sky, T a k y (/, x) 
at a frequency / and position x, consists of several compo- 
nents 

T sky (/, x) = T C mb(/, x) + T gal (/, x) + T ps (/, x) + T m (/, x) , 

(32) 

where Tcmb is the CMB temperature, T ga i is diffuse Galac- 
tic radiation, T ps is emission from (mostly) unresolved extra- 
galactic sources, and Thi is the HI emission we are intending 
to detect. If the telescope and receiver are perfect, the total 
system temperature is given by T sys = T a k y + Tnoise i-e. no 
other contributions (e.g. from ground-spillover) will add to 

Tsys- 

The contribution from each foreground can be decom- 
posed into a smooth component, T(f), and frequency and 
position varying fluctuations, 5T(f, x), about this smooth 
background 

r(/,x)=T(/) + 5T(/,x). (33) 

The foregrounds will be much brighter than the HI sig- 
nal, by several orders of magnitude. At / = 1 GHz, T s k y ~ 
5K, while the HI brightness temperature is ~ 0.1 mK. 
However, our objective is to measure fluctuations of the 
sky brightness temperature as a function of both angular 
scale and redshift. Fluctuations of the HI signal are ex- 
pected to be at the same level as the smooth component 
(STm ~ Thi ~ 0.1 mK), while the fluctuations of the fore- 
grounds will typically be much less (ST <C T). Furthermore, 
we can choose the coldest regions of sky to minimise the 
foreground signal. We will show that regions exist where the 
fluctuations in the foregrounds are in fact a little less than 
two orders of magnitude brighter than the HI signal. Given 
that the foregrounds will be smooth in frequency, this gives 
us confidence that foreground separation algorithms (e.g. 
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(Liu & Tegmark 2011 1) can be employed that will allow a 



detection of BAO. 

Estimates of the foreground brightness temperature 
fluctuations are presented in the next sections for each com- 
ponent and are summarised in Table [2] 



4.2 Galactic synchrotron radiation 

Synchrotron radiation is emitted by electrons spiralling in 
the Galactic magnetic field. The volume continuum emissiv- 
ity is given by the integration of the power radiated by an 
electron P(f; E) at a given energy E and frequency /, and 
the cosmic ray energy spectrum N C (E) 



4tt 



P(f;E)N e (E)dE. 



(34) 



The integration over a distribution of electrons, which itself 
is smooth, results in synchrotron spectra that are smooth, a 
property that we can exploit to separate synchrotron emis- 
sion from hydrogen line emission. 

At frequencies ~ 1 GHz observations indicate that syn- 
chrotron emission is well-approximated by a power-law with 
spectral indices in the range —2.3 to —3.0 (Reich & Re 
~~ 19881). The mean value /3 



ich 



-2.7 and mean dispersion 



ap = 0.12 are measured by Platania et al. (20031 in the 
range 0.4-2.3 GHz. 

A number of large-area surveys are available at frequen- 
cies of ~ 1 GHz, most notably the all-sky 408 MHz map 
of ( jHaslam et~al~1|1982| >. We use the NCSA destriped and 
source-removed version of this marj^Jto estimate the contri- 
bution of synchrotron emission at 1 GHz. The 408 MHz map 
is extrapolated to 1 GHz using a spectral index of /3 = —2.7. 
The map is shown in RA-Declination coordinates in Fig. |13| 
The Galactic plane is evident as well as large-scale features 
that emit over large areas of the high latitude sky. The most 
prominent feature is the North Polar Spur which stretches 
over ~ 100° of the northern sky. The mean temperature at 
high Galactic latitudes (|6| > 30°) is T s k y ~ 2K with r.m.s. 
fluctuations in maps of 1° resolution of <5TG a i ~ 530 mK. The 
r.m.s. fluctuations are dominated by contributions from the 
large-scale features. 

Fortunately, there exist large areas of low synchrotron 
temperature in the north. The declination strip at 8 ~ +40° 
is well known to be low in foregrounds ( jDavies et al. 1(1996] ) 
and is accessible by observing sites in the north. Further- 
more, we are interested in the fluctuations at a given angu- 
lar scale, or multipole value £, rather than the r.m.s. value 
in the map at a given resolution. The r.m.s. estimate from 
the map will be considerably larger due to fluctuations on 
larger angular scales than the beam. 

To estimate the foreground signal at a given lvalue we 
calculated the power spectrum of the 8 = +45° strip af- 
ter masking out the Galactic plane (|6| > 30°). This was 
achieved by using the PolspiceQ software, which calculates 
the power spectrum from the 2-point pixel correlation func- 
tion ( jSzapudi et~aT1|2001| |Chon et all]|2004| ). We apodized 
the correlation function with a Gaussian of width a — 10° to 
reduce ringing in the spectrum at the expense of smoothing 



4 http://lambda.gsfc.nasa.gov/ 

° http: //www2 . iap.fr/users/hivon/software/PolSpice/ . 




Figure 13. In the top plot, we present the entire synchrotron sky 
in celestial (RA/Dec) coordinates, with RA= 0° at the centre and 
increasing to the left. The map has been extrapolated from the all- 
sky 408 MHz ma p ||Haslam et al |1982| | to 1 GHz using a spectral 
index /3 = —2.7 (Pla tania et al.||2003| l. The colour-scale is asinh 
from IK to 8K. The saturated pixels (T > 8K) are shown as 
grey. In the bottom plot, we present three drift scan strips, each 
10° wide, located at declinations of 35° (bottom), 45° (middle) 
and 55° (top). 



the spectrum. A power-law fit was then made to estimate 
the foreground signal at I = 200 (~ 1°). The result is given 
in Table H 

For a 10° wide strip centred at 8 = 45° we find that the 
Thai = 1.7K with fluctuations of STq^i = 67 mK. Utilising 
these cool regions will be important for detecting the faint 
HI signal. The exact choice of strip will be a critical factor to 
be decided before construction of a fixed elevation telescope 
while the Galactic mask to be applied can be chosen during 
the data analysis step. 

4.3 Free-free radiation 

Free-free radiation (thermal bremsstrahlung) is emitted by 
electrons interacting with ions. At radio frequencies, this 
comes from warm ionized gas with typical temperatures of 
T e ~ 10 4 K. Away from the Galactic plane, the radiation is 
optically thin and its brightness temperature is given by 



90 mK (&) 



' f 

> GH: 



~ 2 1 / EM \ 
/ \ cm - e pc J 



(35) 



where EM is the emission measure representing the in- 
tegral of the electron density squared along the line-of- 
sight (EM= J nidi). The free-free spectrum is therefore a 
well-defined power-law with a temperature spectral index 
/3 — —2.1 which acts to flatten the spectral index of the total 
continuum where it has a brightness temperature compara- 
ble to that of the synchrotron emission. 

The EM can be estimated using Ha measurements, 
which can then be converted to radio brightness tempera- 
ture assuming T e is known and the dust absorption is small 
( Dickinson et al.||2003| ). For \b\ > 30°, where T c » 8000K, 
EM ~ 1.3 cm -0 pc at 1 GHz, this corresponds to Tg ~ 5 mK. 
For the 10° wide strip at 8 — +45° the free-free brightness 
temperature fluctuations are 8Tg = 0.25 mK at I = 200. The 
free-free component is therefore considerably weaker than 
the synchrotron component across most of the sky, but will 
act to (slightly) flatten the spectrum at higher frequencies. 
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Table 2. Summary of foregrounds for HI intensity mapping at 1 GHz for an angular scale of < 
10°-wide strip at declination S = +45° and for Galactic latitudes |6| > 30°. 



200). The estimates are for a 



Foreground 


T 


ST 




[mK] 


[mK] 


Synchrotron 


1700 


67 


Free- free 


5.0 


0.25 


Radio sources (Poisson) 




5.5 


Radio sources (clustered) 




47.6 


Extragalactic sources (total) 


205 


48 


CMB 


2726 


0.07 


Thermal dust 




~ 2 x 10" 


Spinning dust 




~ 2 x 10" 


RRL 


0.05 


3 x 10" 3 


Total foregrounds 


~ 4600 


~ 82 


HI 


~ 0.1 


~ 0.1 



Notes 



Power-law spectrum with j3 RJ —2.7. 
Power-law spectrum with j3 ~ —2.1. 
Assuming removal of sources at S > 10 mjy. 
Assuming removal of sources at 5" > lOmJy. 
Combination of Poisson and clustered radio sources. 
Black-bo dy spectrum, (0 = 0). 



6 Model of Finkbeiner et al. 



119991 



|Davies et al| ( |2006| and CJNM model of 
Hydrogen RRLs with An = 1. 



Draine & Lazarian 



(1998) 



Total contribution assuming the components are uncorrelated. 
Cosmological HI signal we are intending to detect. 



4.4 Extragalactic foregrounds 

Extragalactic radio sources are an inhomogeneous mix of 
radio galaxies, quasars and other objects. These contribute 
an unresolved foreground^] The sources are often split into 
two populations based on their flux density spectral index, 
a (S oc v a ): steep spectrum sources with a < —0.5 and 
flat spectrum sources with a > —0.5. At low frequencies 
(< 1GHz), steep spectrum sources with a as —0.75 will 
dominate ( |Kellermann et al.|1968[ |. 

The contribution of point sources, T ps , can be calcu- 
lated from the differential source count, dN/dS, represent- 
ing the number of sources per steradian, N, per unit flux, 
S. A number of compilations of source counts are available. 
We chose to use data collected from continuum surveys at 
1.4 GHz data between 1985 and 2009, summarised in Ta- 
ble [3] Fig. [14] shows the Euclidean normalized differential 
source count for our compilation. We have fitted a 5th order 
polynomial to these data of the form 



log 



S 2 - 5 dN/dS\ 



No 



log 



So 



(36) 



where iVo = Uy sr and So = Uy are normalizing con- 
stants. A least-squares fit gives the best-fitting coefficients: 
a = 2.593, oi = 0.093, a 2 = -0.0004, a 3 = 0.249, 
ci4 = 0.090 and as = 0.009. These values are compatible 



with those of Vernstrom et al. (2011 I 



Assuming we can subtract all sources with flux density 
S > S m ax, we can calculate the mean temperature 



dB 
dT 



(37) 



where dB/dT = 2fc/ 2 /c 2 , c is the speed of light and k is 
the Boltzmann constant. We find T pB = 97 mK at 1.4 GHz 
for S'max = oo. We would expect to be able to subtract 
the brightest radio sources down to S'max = lOmJy using 
the NR AO VLA Sky-Survey (NVSS) which is complete to 
3.4 mjy ( Condon et al.|1998[ ) although there will be a resid- 
ual contribution due to the variability of radio sources. Sub- 



6 For convenience we refer to extragalactic sources as a fore- 
ground though both lower redshift (foreground) and higher red- 
shift (background) sources will contaminate the HI signal. 



Table 3. References used for source surveys at 
for compilation of source counts. 



1.4 GHz used 



1! cfcrcncc 



jBondi et aT||2003| 
[Ciliegi et^L](l999) 
Fomalont eTal!|] 2006 
^ruppioni et al . ( 1999 
Hopkins et al. ( 1999 
|lbar et aT^2009| 
%THondoirri[l985 
fc Morrison| p OUf 
Richards| l|2000| 
Seymour et al. ( 2008 ) 



n 

:»6) 

99) 
9) 



White et al. 



(2008) 



ftCiliegi et oL (1999) 
O Fomalont et al. (2006) 
iGruppioni M tf. (1989) 
□ Hopkins et ol. (2003) 
v icar et al. (2009) 
O Mitchell Ji Condon (1985) 
*Oen ft Morrison (2008) 
V Richards (2000) 
• Seymour et al. (2008) 




*** 



5 (mjy) 

Figure 14. The compilation of differential radio sources counts 
at 1.4 GHz which we have used. The solid line is the best-fitting 
5th order polynomial. 



tracting sources stronger than 10 mjy reduces the mean 
temperature of T ps = 83 mK at 1.4 GHz. Although this is 
a useful thing to do, we clearly see that the bulk of the 
background temperature arising from extragalactic sources 
is contributed by the weakest ones and this is precisely the 
point where the uncertainty in the counts is greatest. Ex- 
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trapolating the course counts to 1 GHz with the mean tem- 
perature spectral index f3 = — 2.7 gives T ps = 205 mK. 

Before we continue, we note that there is evidence of 
an additional component to the sky brightness as measured 
by the ARCADE experiment at 3 GHz ( |Fixsen et al.|2011 |. 
They measured an excess component of radiation with a 
brightness temperature of ~ 1 K at 1 GHz and with a spec- 
tral index /3 ~ —2.6. Given the extreme difficulty in mak- 
ing absolute sky brightness measurements further observa- 
tions, perhaps at lower frequency where the claimed excess 
is stronger, are required to firmly establish its reality. How- 
ever, if it is real, this could require either a re-interpretation 
of Galactic radiation contributions or be evidence for a new 
population of extragalactic sources ( Seiffert et al.|2011| > - cf 
our point concerning the sensitivity of the computation of 
T ps to the lower cut-off. For our initial analysis presented 
here, we do not consider this result further although we 
note that this could increase our foreground estimates sig- 
nificantly. 

There are two contributions to ST ps . The first is due 
to randomly (Poisson) distributed sources. However, if the 
sources have a non-trivial two-point correlation function, 
then there is a contribution due to clustering. We quantify 
the spatial fluctuations in ST ps in terms of the angular power 
spectrum, C t = cf oisson + cf ustCTins , which is the Legendre 
transform of the two point correlation function C(9). 

Poisson distributed sources have a white power spec- 
trum given by 



C 



dB\- 2 
dT ) 



C 2 dN 

s ds dS 



(38) 



In terms of multipoles, I (where 9 ~ 180°/£), the magnitude 
of the temperature fluctuations at any given angular scale is 
given by ST = \J i(l + l)Ci/2ir. Using our source count fit 
above to integrate up to S'max = lOmJy and extrapolating 
to 1 GHz leads to fluctuations of 5T pB = 5.5 mK at I = 200 
(ldeg) 

The clustered component of radio sources acts to in- 
crease pixel-pixel correlations. The power spectrum due the 
clustered sources can be simply estimated as cf UBtCT = 
w eTp S , where we is the Legendre transform of the angular 



correlation function, w{9) (Scott fc White||T9 99 ) . The clus- 
tering of radio sources at low flux densities (< lOmJy) is 
not well-known. To make an estimate, we use w(9) mea- 
sured from NVSS, w hich can be approxim ated as w(9) m 
(1.0±0.2) x 1O" 3 -0 - 8 (Overzieretal.|2003| ). Legendre trans- 
forming this gives wt ~ 4.8 x 10~ i £~ L ' 1 . Although measured 
directly for sources at a higher flux density than those rele- 
vant here, the overall amplitude and power law are compat- 
ible with those expected if the galaxies follow clustering of 
cold dark matter. 

The contribution to ST ps at 1 GHz on angular scales 
~ 1° from clustered sources, assuming a cut-off flux density 
of lOmJy, is 8T ps w 47.6 mK. Plotting both the Poisson and 
clustering contributions to ST ps , it can be seen in Fig. |I5| that 
the contribution from the clustering component appears to 
be dominant over the Poisson fluctuations for multipoles 
£ < 500 and is approximately an order of magnitude above 
the fluctuations in the HI signal. 

The total contribution to ST ps is then the quadrature 
sum of the Poisson and clustered components (they are un- 
corrected), i.e. Cf = cf oisson + cf UBtclins which leads 



S„„=10 mjy 



!-o:ssnn 
Clustering 




Figure 15. The Poisson and clustering contribution to the sky 
brightness temperature fluctuations for a flux cut-off value of 
S'max = lOmJy. The brightness temperature fluctuations due to 
the 21cm signal is expected to be < 300 /^K at redshifts z S3 0.5 
and is indicated by the hatched region. 



to ST B . 



\/^^Poiaaon + ^cluster 



48 mK. We note that 



there is considerable uncertainty in this estimate, particu- 
larly with the clustered component. The extrapolation of 
our source count into the low flux density (<C 0.1 mjy) 
regime may be unreliable (e.g. because the ARCADE re- 
sults may imply an unexpectedly large population of very 
weak sources) and can only be firmly established by ultra- 



deep radio surveys with upcoming telescopes (e.g. |Norris et 



(20111) 



4.5 Other foregrounds 

A number of other foregrounds could also contribute to T s k y - 
Other diffuse continuum emission mechanisms include ther- 
mal dust and spinning dust emissions. Both of these are 
expected to be negligible at v ~ 1 GHz. For example, ex- 
trapolation of the thermal dust model of |Finkbeiner et al.| 
( 1999[ ), at 1° resolution predicts 6T ~ 2 x 10" e mK for 
\b\ > 30°. Similarly, extrapolating the anomalous microwave 



emission observed at 23 GHz by WMAP ( |Davies et al.|2006 \ 
and assuming it is all due to spinning dust emission from the 
cold neutral medium (CNM; Draine fc Lazarian|1998 1, gives 



ST ~ 2 x 10- 5 mK and can therefore be safely neglected. 

There are a number of spectral lines which emit in the 
radio band. The most prevalent in the frequency range of in- 
terest are radio recombination lines (RRLs), primarily from 
hydrogen. The frequencies of these lines are given by the Ry- 
dberg equation, for example, H187q is at a rest frequency 
of 997.63 MHz. The intensity of these lines can be estimated 
from the free-free continuum temperature, Tg, using the fol- 
lowing equation ( Rohlfs fc Wilson|2004 1 

jndv 

T H AV 



6.985x10" 



a(T e 



l + n(H) /T c 
n(He) \K 



GHz 



(39) 

where a(T e , v) is a correction factor, which at low frequencies 
is very close to 1 ( Mezger fc Henderson|1967 1 . For a typical 
linewidth of 25kms~ , assuming T e = 8000 K and a 8% he- 
lium to hydrogen ratio, the line temperature will be ~ 1 % 
of the free-free continuum temperature at 1 GHz. Given the 
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estimates for Galactic foreground emission from Section |4.3| 
we expect f L ~ 0.05 mK and ST L ~ 2.5 x 10" 3 mK. There- 
fore, RRLs are not expected to be an important foreground 
either. Nevertheless, it would be wise to omit frequency 
channels that are known to contain the brightest (An = 1) 
hydrogen lines, which occur at intervals of approximately 
A/ S3 20 MHz. We note, however, that these recombination 
lines contain much astrophysically interesting information 
about the Galaxy and therefore are likely to provide signif- 
icant spin-off science from any radio BAO project. 



5 THE BINGO CONCEPT 

The discussion presented in Section [3] suggests that in order 
to get a good detection of BAO at low redshifts (case A from 
the optimisation study) a suitable single telescope should 
satisfy the following design constraints: 

• an angular resolution defined by #fwhm ~ 40 arcmin 
corresponding to an under-illuminated primary reflector 
with an edge taper < — 20 dB; 

• operate in the frequency range above 960 MHz with a 
bandwidth that is as wide as possible - exactly how wide 
and exactly where the band is positioned will ultimately be 
determined by those parts of the spectrum relatively free of 
radio frequency interference (RFI); 

• frequency resolution should be ~1 MHz or better; 

• frequency baseline should be as free as possible from 



Scanning direction y 



instrumental ripples, particularly standing waves (Rohlfs & 
|Wilson|[2TJ04t ; 

• the receivers should be exceptionally stable in order to 
allow thermal-noise-limited integrations of a total duration 
long enough to reach a surface brightness limit ~ 100/iK per 
pixel; 

• the sky coverage as large as possible, but at the very 
least >2000 deg 2 with a minimum size in the smallest di- 
mension of 10 deg; 

• the number of feeds should be tjf > 50 in order to 
achieve a sufficiently large signal-to-noise ratio; 

• sidelobe levels as low as possible so as to minimize the 
pickup of RFI and emission from strong sources; 

• beam ellipticity < 0.1 in order to allow map-making 
and power spectrum analysis routines to work efficiently; 

For reasons of stability, simplicity, and cost, a telescope 
with no moving parts would be very desirable. We suggest 
that the following conceptual design could deliver what is 
required. We first concentrate on the telescope and then 
on the receiver design. We have named the overall project 
BINGO which stands for BAO from Integrated Neutral Gas 
Observations. 



5.1 The telescope and focal plane 

Operation at frequencies between 960 MHz and 1260 MHz 
is a practical possibility because it avoids strong RFI from 
mobile phone down-links which occupy the band up to 960 
MHz. In order to achieve the desired resolution (better than 
40 arcmin at 1 GHz), an under-illuminated primary reflector 
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Figure 16. Focal plane of staggered feeds. A total of 60 pixels 
laid out in a rectangle (dashed black line) of 16 mx 15 m. This will 
allow for a FOV of about 10x9 deg 2 . Axes x and y are shown as a 
reference to the positions mentioned in Table [4] The pixel shown 
in red is the off-axis pixel corresponding to the beam profiles in 
Fig.[T7| 



of about 40 m projected aperture diameter, D, is required^] 
If the telescope is to have no moving parts, i.e. fixed pointing 
relative to the ground, then a wide instantaneous field-of- 
view (FOV) with multiple feeds is essential so that as the 
Earth rotates a broad strip of sky can be mapped during the 
course of a day. To reach the target sensitivity, a minimum 
number of 50 pixels is required in the focal plane, each of 
them reaching the following performance: 

• Forward gain loss in comparison to central pixel < 1 dB; 

• Cross-polarisation better than — 30 dB; 

• Beam ellipticity lower than about 0.1. 

We propose to adopt an offset parabolic design in order 
to avoid blockage, minimise diffraction and scattering from 
the struts, and also to avoid standing waves that might be 
created between the dish and the feeds. A wide FOV with 
the required optical performance mentioned above can only 
be achieved with a long focal length. We found that a fo- 
cal length F — 90 m, giving a focal ratio F/D = 2.3, much 
larger than more conventional radio telescopes (having typ- 
ically a F/D ratio of about 0.4), is suitable for our purpose. 
This has the drawback of requiring the use of feeds with a 
small full-width-half-maximum, therefore with a large aper- 
ture diameter. For our telescope design, this means feeds 
with an aperture diameter of at least 2 m producing an edge 
taper telescope illumination of —20 dB in order to reduce the 
spillover. We are targeting a spillover of less than 2% in to- 
tal power as the fraction of the horn beam not intercepting 
the telescope will see radiation from the ground therefore 
increasing the system temperature T aya . 

The design adopted for the focal plane is given in 
Fig. [16] A line of eight feeds allows for coverage of a 10 deg 

7 The equivalent fully illuminated aperture would be D S3 25 m 
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Table 4. Optical performance calculated using GRASP assuming 
ideal telescope and Gaussian feeds with -20 dB edge taper. Results 
are given for a pixel that would be located at the centre of the 
focal plane and for the corner edge pixel (red in Fig. |16| at x = 
-8m;j = —7.5 m). Peak cross-polarisation is given respectively 
to the co-polarisation forward gain. 



Pixel 


Forward 


Ellipticity 


Peak 


FWHM 


Position 


Gain 




X-Pol 






dB 




dB 


arcmin 


Centre 


49.8 


io~ 3 


-40 


38 


Edge 


49.6 


0.07 


-35 


39 



Table 5. Summary of proposed BINGO telescope parameters. 



Reflector diameter 


40 m 


Resolution (at A = 30 cm) 


40 arcmin 


Number of feeds 


60 


Instantaneous field of view 


10 degx9 deg 


Frequency range 


960 MHz to 1260 MHz 


Number of frequency channels 


> 300 



strip. Along the scanning direction produced by the Earth's 
rotation, other lines of feeds are staggered to allow for a 
proper sampling on the sky. With this configuration, 60 pix- 
els, 10 more than the requirement, will produce a 10 x 9 deg 2 
FOV. Assuming a feed aperture diameter of about 2 m, the 
resulting focal plane will be comprised within a rectangle of 
16 m x 15 m. 

Optical simulations of such a concept have been per- 
formed using Physical Optics (GRASP from TICRA) in 
order to take into account the diffraction effects. Results 
from this model, assuming ideal components (perfect tele- 
scope and Gaussian feedhorns) and a perfect alignment, are 
given in Table|4] The performance is well within the require- 
ments even for the most off-axis pixel (represented in red in 



Fig. 161. Fig. 1 17| represents the modelled contour plot of the 
beams (co- and cross-polarisation beams) in the U-V plane 
that can be expected for the pixel located at the corner edge 
of the focal plane. Taking into account degradations due to 
using real components and with typical alignment errors, we 
can till expect to meet the requirements. 

The design concept is illustrated in Fig |18| and the key 
parameters listed in Table [5] Rather than a high tower to 
support the focal arrangement, the plan would be to identify 
a suitable site in which the dish structure can be built below 
a cliff (or near vertical slope) of height 90 m so that the feed 
arrangement can be fixed on a boom near the top of the cliff. 



5.2 Design for a receiver module 

A high degree of stability is an essential design criterion 
both for telescope and receivers. Such demands are not 
unique since instruments used to measure the structure in 
the CMB have been built which meet similar stability cri- 
teria. In CMB experiments correlation receivers (or inter- 
ferometers) are generally the design of choice (for example, 
WMAP and Planck). They can reduce the effect of gain fluc- 
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Figure 17. U-V maps of beam intensity for the most off-axis 
pixel corresponding to the focal plane position x = —8m and 
y = —7.5m (represented in red in Fig |16[ l. Top: co-polarisation; 
bottom: cross-polarisation. Simulations performed using ideal 
telescope and horn. Colour scale: the gain in dB relative to max- 
imum. 



tuations in the amplifiers (1// noise) by up to three orders 
of magnitude ( Mennella et al.|2010 |. 

We plan to adopt the same approach with each receiver 
module producing a "difference spectrum" between the de- 
sired region of sky and a reference signal. The choice of the 
reference is important. Ideally it should have the same spec- 
trum as the sky and be close to the same brightness tem- 
perature. It must also be very stable. An obvious choice is 
another part of the sky. We propose that there should be 
a reference antenna for each module with a beamwidth of 
around 20° which would have a beam area sufficiently large 
to average out all but the lowest I modes (see Fig. [2|. It 
should be fixed and we suggest that it should point directly 
to one of the celestial poles. In this way all measurements at 
all times of day, and all the year round, would then be com- 
pared to the same region of sky. Since the sky emission has 
significant linear polarization (~ 10 % or more; Wolleben 
et al. (2006 1), measuring with circularly polarized feeds is 



necessary. Similar considerations will be important for fore- 
ground emission removal. A schematic diagram for a single 
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Figure 18. Schematic of the proposed design of the BINGO telescope. There will be an under-illuminated ~ 40 m static parabolic 
reflector at the bottom of a cliff which is around ~ 90 m high. At the top of the cliff will be placed a boom on which the receiver system 
of 60 feed-horns will be mounted. 
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Figure 19. Block diagram of the receiver chain. The reference beam will point towards one of the Celestial Poles. The hybrids will be 
waveguide magic tees. After the second hybrid the signal will be down-converted using mixers with a common local oscillator. There will 
be further amplification before digitization. 



polarization channel of the proposed receiver module design 
is shown in Fig. [19] 



5.3 Receiver simulations 

An ideal correlation receiver - one with perfectly matched 
system temperatures and no imbalance in the power split- 
ting within each hybrid - can effectively remove all 1// noise 
from its signal output. If A = A sys + 5 A and B = B sys + SB 
are the total voltages entering the receiver chain and A sya oc 



/T^andB sys 



oc \fT^ s are the system temperature con- 
tributions and S A and SB are the signals we want to detect, 
then the difference detected power at the end of the chain 
is proportional to the difference, 8 A 2 — SB 2 , between the 
science (A) and reference beam (B). The spectrum of such 
a signal exhibits no excess power at low frequencies (1// 

and integration can con- 



noise) as illustrated in Fig. 20(a 



tinue indefinitely with the white noise falling as the square 
root of integration time. 

However, real life receivers are imperfect. Imbalances 
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between system temperature and power splitting both limit 
the reduction of the 1/f noise resulting in a break in the 
power spectrum quantified in terms of a knee frequency 
1 Seiffert et al. |2002 \ separating white noise-dominated and 
pink noise-dominated parts of the spectrum. Consequently 
there will be a limit to the viable integration time as illus- 



trated in Fig. 20(b) It is useful to quantify by how much a 



correlation receiver is limited by these imperfections, and to 
do so wc simulate their effect in a realistic receiver model. 
We assess the maximum permissible imperfection amplitude 
to ensure the knee frequency remains below /knee ~ 1 mHz, 
a value corresponding ~ 20 minutes integration time - suffi- 
cient for the largest HI structures we are interested in detect- 
ing, ~ a few degrees, to drift through the telescope beam. 

We include Gaussian (white) and 1/f noise, generating 
the former using a simple random number procedure and 
the latter with a more complex auto-regressive approach - 
essentially filtering a white noise time stream into one with 
correct correlation properties. The knee frequency of com- 
ponent amplifiers depends on both the physical choice of 
amplifier and the channel bandwidth. Though the actual 
value will need to be determined empirically, we will assume 
in the current simulations that the individual amplifier knee 
frequencies are /^™^ ~ 1 Hz for a channel bandwidth of 
1MHz. 



5.3.1 Unbalanced T SJ 

In an ideal receiver th 
be perfectly matched. This is impossible to achieve at all 
times due to the angular variations in the Galactic emission 
and here we investigate the effect on the knee frequency due 
to small imbalances between the temperatures. 

Modelling the entire correlation receiver and simulat- 
ing noise spectra for the output difference signal in the 
presence of small temperature imbalances, we find that the 
knee frequency of the correlation receiver can be approxi- 
mated by /knee = 1.2Hz((A sy3 — B sys )/A sys ) 2 and there- 
fore to ensure BINGO 's knee frequency remains continu- 
ally below the desired /knee ~ 1 mHz, our simulations show 
that (A sys — B sys )/ A sy3 needs to be < 0.03. In practice the 
balance can be maintained electronically so this should not 
be an important limitation, provided the unbalance behaves 
predictably as a function of time. The main source of imbal- 
ance in the system temperatures will arise from the changing 
sky background as the telescope scans across the sky. Typi- 
cal variations in those parts of the sky of interest well away 
from the Galactic plane will be a few Kelvin and will be re- 
peatable from day to day. Importantly the sky emission will 
be spectrally smooth and virtually identical in all frequency 
channels. This means that once the frequency-averaged vari- 
ation of sky brightness has been measured over the course 
of a siderial day, a multiplicative correction, the r-factor de- 
fined by dSeiffert et al. |2002| >, can be applied to all channels 
to bring the A and B outputs into balance before differenc- 
ing, thus ensuring the (A sys 



B sy3 )/A sys < 0.03 is met. 




5.3.2 Power Split Imbalance 

We now consider the case of hybrids in which the power 
splitting is imperfect. To do so, we introduce the constant 



coefficient e, representing an imbalance between the power 
emerging from the two output ports of the hybrids. The 
signals emerging from such a hybrid are then 



(40) 



where A, B are the two input signals. 

Here we assume the power splitting imbalances in the 
first and second hybrids are equal, and simulate noise spec- 
tra for different values of e. Fitting for the knee frequency 
we find /knee = 4e 2 Hz. Again, in order to maintain a knee 
frequency of 1 mHz, this relationship implies that the power 
splitting imbalance should satisfy e < 0.015. This require- 
ment can be relaxed if a 180° phase switch is used within the 
hybrid loop and, since waveguide magic-tees typically have a 
power split imbalance of e = 0.02, the required performance 
should be achievable. 



5.4 Practicalities 

Our philosophy is to choose simplicity while keeping in mind 
good engineering practice and, wherever possible, to adopt 
conventional solutions to problems rather than try to invent 
new ones. 

To reduce the impact of RFI the telescope would need 
to be in a remote site and therefore an important additional 
consideration is that everything should be easy to maintain. 
The design, construction and mounting of the feeds is per- 
haps the most challenging aspect of the overall concept. As 
stated above, low sidelobe levels are important in order to 
mitigate the effects of RFI (even in the remotest part of 
the World) and ease the separation of foreground emission 
from the desired neutral hydrogen signal. If a conventional 
horn design is adopted, these horns will have to be extremely 
large; to illuminate a 30 m dish from a focus ~ 75 m away 
will require an aperture diameter for each feed of 1.5 to 2 m 
and a length of ~ 8 m. Such feeds are practical, in the sense 
that they would be scaled up versions of horns built to ex- 
isting feed designs, but could prove expensive to build and 
unwieldy. Alternatives like using small secondary mirrors for 
each receiver module are under consideration. 

An integral part of the feed system is the extraction 
of polarized signals which we will do by extracting the two 
opposite circular polarizations. Minimizing any Ohmic loss 
before the first amplifiers is essential if good thermal noise 
performance is to be achieved. After the feed the next ele- 
ment is the hybrid. Due to their low loss and good power 
splitting properties we would use waveguide magic-tees. 

We would use non-cryogenic low noise amplifiers be- 
cause at this frequency the noise advantages of cooling to 
~ 20 K are outweighed by the low cost, the low power 
consumption and the low maintenance demands of a non- 
cryogenic approach. It would also be desirable to integrate 
the amplifiers with the magic-tee to avoid the necessity to 
use potentially lossy cables. Some cooling, perhaps using 
Peltier coolers, could be used because it would provide both 
temperature stabilization for the amplifiers and a potential 
improvement in noise performance. 

An important choice is when to digitize the signal. In 
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Figure 20. Simulated noise power spectra for a realistic BINGO receiver which uses amplifiers with /knee = 1 Hz and a 30 min integration. 
Here (a) represents a perfectly balanced receiver; and (b) a receiver with imbalances in both system temperature, (A Bya — B Bys )/A syB = 0.2 
and power splitting, e = 0.4. These values are significantly larger than expected in a realistic situation and have been chosen to illustrate 
the point. The upturn of the spectrum in (b) indicates the presence of 1// noise. 



conventional correlation receivers the "correlation" is per- 
formed in a second hybrid. Our preferred solution is to down- 
convert the radio frequency signal with a common local os- 
cillator and then to digitize and perform the function of 
the second hybrid in software using FPGAs to multiply the 
signals and produce averaged spectra which can then be dif- 
ferenced. One advantage of such a system is that the levels 
in the two arms of the receiver can be easily adjusted so that 
they are the same and this then operates the receiver at its 
optimum for gain fluctuation reduction. 



significantly affect our ability to extract the HI signal due to 
the fact that its spectral profile is expected to be relatively 
smooth. 

We have put forward a design concept, which we have 
named BINGO, to perform the lower of the two redshift 
surveys studied. This will use a large focal length tele- 
scope, F/D ~ 2.3, to mount an array of 60 feedhorns ap- 
proximately 90 m above a static 40 m. It will use pseudo- 
correlation receivers to reduce 1/f noise due to gain fluctu- 
ations. We are currently searching for a site and the funding 
needed to bring this idea to fruition. 



6 CONCLUSIONS 

Intensity mapping is an innovative idea which has been sug- 
gested to allow large-scale structure surveys using the 21cm 
line of neutral hydrogen using instruments with relatively 
small collecting area compared to the SKA. Such concepts 
would naturally be much cheaper. We have presented calcu- 
lations that illustrate the characteristics of the signal and it 
is shown that it should have an r.m.s. of ~ 100 fiK on an- 
gular scales ~ 1 deg and bandwidths of 1 — 10 MHz. At very 
low redshifts there will be a shot noise effect due to it being 
dominated by individual galaxies, at medium redshifts it is 
dominated by the effects of HI in galaxy clusters and at high 
redshifts becomes an unresolved background. 

We have performed an optimisation analysis which sug- 
gests that a significant measurement of the acoustic scale can 
be made between z = 0.13 — 0.48 using a single dish tele- 
scope of diameter 40 m with a drift scan survey of a strip 
with width 10 deg covering a total of ~ 2000 deg 2 . In or- 
der to do this at higher redshift (specifically, the case of 
z — 0.58 — 1.47) would require a much larger dish, ~ 140 m, 
and an area of 500 deg 2 which would probably be best done 
using an interferometer in order to allow non-drift scan ob- 
servations using such a long baseline instrument. We have 
presented a detailed assessment of the likely foregrounds and 
conclude that there will be continuum fluctuations with an 
r.m.s of ~ 100 mK on degree scale, but that this should not 



ACKNOWLEDGEMENTS 

We thank Peter Dewdney for suggesting the use of a cliff 
to mount the BINGO feeds and Xiang-Ping Wu for sug- 
gesting using one of the Celestial Poles as a reference. We 
have had numerous useful discussions with Peter Wilkinson, 
Shude Mao and Richard Davis. CD acknowledges an STFC 
Advanced Fellowship and an EC IRG grant under the FP7. 
AP is supported by the project GLENCO, funded under 
the FP7, Ideas, Grant Agreement n. 259349 and also re- 
ceived support from the STFC(UK). We thank Bob Watson 
for help in creating Fig. |18| 



REFERENCES 

Abbott, T. et al., 2005, arXiv:astro-ph/05 10346 
Abdalla, F. B. and Rawlings, S., 2005, MNRAS, 360, 27 
Abdalla, F. B., Blake, C. and Rawlings, S., 2010, 401, 734 
Afshordi, N. and Loverde, M., 2008, Phys. Rev. D78, 
123506 

Anderson L. et al., 2012, arXiv: 1203 . 6594 
Battye, R. A., Davies, R. D. & Weller, J., 2004, MNRAS, 
355, 1339 

Beutler, F., Blake, C, Colless, M., et al. 2011, MNRAS, 
416, 3017 



HI intensity mapping : a single dish approach 19 



Blake, C. A., Davis, T., Poole, G., et al. 2011, MNRAS 
415, 2892 

C. Blake and K. Glazebrook, Astrophys. J. 594 2003) 665 
K. Glazebrook & C. Blake, 2005, Astrophys. J. 631 1 
Bondi, M., Ciliegi, P., Zamorani, G., et al. 2003, A&A, 403, 

857 

Bonvin, C, Durrer, R., 2011, Phys. Rev. D84 063505 
M. Bowden et al, Mon. Not. Roy. Astron. Soc. 349 (2004) 
321 

E. F. Bunn and M. J. . White, Astrophys. J. 480 (1997) 6 
S. M. Carroll, W. H. Press and E. L. Turner, Ann. Rev. 

Astron. Astrophys. 30 (1992) 499. 
Chang, T. -C, Pen, U. -L., Peterson, J. B. et al. 2008, 

Phys. Rev. Lett., 100, 1303 
Challinor, A., Lewis, A., 2011, Phys. Rev. D84 043516 
Chang, T. -C, Pen, U. -L., Bandura, K. et al. 2010, Nature, 

466, 463 

Chon, G., Challinor, A., Prunet, S., Hivon, E., & Szapudi, 
I. 2004, MNRAS, 350, 914 

Ciliegi, P., McMahon, R. G., Miley, G., et al. 1999, MN- 
RAS, 302, 222 

Clifton, T., Ferreira, P. G., Padilla, C, et al. 2011, 
arXiv: 1106.2476 

Cole, S., Percival, W. J., Peacock, J. A., et al. 2005, MN- 
RAS, 362, 505 

Condon, J. J., Cotton, W. D., Greiscn, E. W., et al. 1998, 
AJ, 115, 1693 

Copeland, E. J., Sami, M., Tsujikawa, S., 2006, IJMPD, 
15, 1753 

Davies, R. D., Watson, R. A., & Gutierrez, C. M. 1996, 

MNRAS, 278, 925 
Davies, R. D., Dickinson, C, Banday, A. J., et al. 2006, 

MNRAS, 370, 1125 
Dickinson, C, Davies, R. D., & Davis, R. J. 2003, MNRAS, 

341, 369 

Draine, B. T., & Lazarian, A. 1998, ApJ, 508, 157 
Duffy, A. R., Battye, R. A., Davies, R. D., et al. 2008, 

MNRAS, 383, 150 
Duffy, A. R., Kay, S. T., Battye, R. A., Booth, C. M., Dalla 

Vecchia, C, Schaye, J., 2012, 420, 2799 

D. J. Eisenstein and W. Hu, Astrophys. J. 496 (1998) 605 
Eisenstein, D. J., Hu, W., Tegmark, M., 1998, ApJ, 504, 

L57 

Eisenstein, D. J., Annis, J., Gunn, J. E., et al. 2001, As- 
trophys. J., 122, 2267 

Eisenstein, D., J., 2003, in Dey, A. (Ed), Wide-field Multi 
Object Spectroscopy, ASP Conference Series 

Eisenstein, D. J., Zehavi, I., Hogg, D. W., et al. 2005, ApJ, 
633, 2 

Feldman, H., Kaiser, N., Peacock, J., 1995 AJ 426, 23 
Finkbeiner, D. P., Davis, M., & Schlegel, D. J. 1999, ApJ, 
524, 867 

Fixsen, D. J., Kogut, A., Levin, S., et al. 2011, ApJ, 734, 5 
Fomalont, E. B., Kellermann, K. I., Cowie, L. L., et al. 

2006, ApJS, 167, 103 
Furlanctto, S. R., Oh, S. P., Briggs, F. H., 2006, Phys. 

Rep., 433, 181 

Giardino, G., Banday, A. J., Gorski, K. M., et al. 2002, 
A&A, 387, 82 

Gruppioni, C, Ciliegi, P., Rowan-Robinson, M., et al. 1999, 
MNRAS, 305, 297 
Haslam, C. G. T., Salter, C. J., Stoffel, H, & Wilson, W. E. 



1982, A&ASS, 47, 1 
Hawken, A. J., Abdalla, F. B., Hiitsi, G., Lahav, O., 2011, 

arXiv: 1111.2544 
Hopkins, A., Afonso, J., Cram, L., & Mobasher, B. 1999, 

ApJL, 519, L59 
Ibar, E., Ivison, R. J., Biggs, A. D., et al. 2009, MNRAS, 

397, 281 

Jarosik, N. et al. 2011, ApJS, 192, 14 

Kellermann, K.I., Paulint-Toth, U.K., Davis. M.M., 1968, 

ApL, 2, 105 
Knox, L. 1995, Phys. Rev. D, 52, 4307 
Komatsu, E. et al. 2011, ApJS, 192, 18 
Laureijs, R. et al. 2011, arXiv: 1110 . 3193 
Limber, D.N., 1953, ApJ, 117, 134 

Liu, A., & Tegmark, M. 2011, Phys. Rev. D, 83, 103006 
Loeb, A. & Wyithe, J.S.B., 2008, Phys. Rev. Lett, 100, 
161301 

Martin, A.M., Giovanelli, R., Haynes, M.P., Guzzo, L., 

2012, ApJ, 750, 38 
Mennella, A., Bersanelli, M., Butler, R. C, et al. 2010, 

A&A, 520, A5 

Meyer, M. J., Zwaan, M. A., Webster, L., et al. 2004, 
MNRAS, 350, 1210 

Mezger, P. G. & Henderson A. P., 1967, ApJ, 147, 471 

Mitchell, K. J., & Condon, J. J. 1985, AJ, 90, 1957 

R. P. Norris, A. M. Hopkins, J. Afonso, S. Brown, J. J. 
Condon, L. Dunne, I. Feain, R. Hollow, M. Jarvis, M. 
Johnston-Hollitt, E. Lenc, E. Middelberg, P. Padovani, 
and other 44 authors, 2011, PASA, 28, 215 

Obreschkow, D., Croton, D., Lucia, G., Khochfar, S., Rawl- 
ings, S., 2009, ApJ, 698, 1467 

Obreschkow, D., Klockner, H. -R., Heywood, I., Levrier, 
F., Rawlings, S., 2009, ApJ, 703, 1890 

Overzier, R. A., Rottgering, H. J. A., Rengelink, R. B., & 
Wilman, R. J. 2003, A&A, 405, 53 

Owen, F. N., & Morrison, G. E. 2008, AJ, 136, 1889 

Pen, U.-L., Staveley-Smith, L., Peterson, J., Chang T.-C. 
2009, MNRAS, 394, L6 

Percival, W. J., Baugh, C. M., Bland-Hawthorn, J., et al. 
2001, MNRAS, 327, 1297 

Percival, W. J., et al. 2002, MNRAS, 337, 1068 

Percival, W. J., Reid, B. A., Eisenstein, D. J., et al. 2010, 
MNRAS, 401,2148 

Perlmutter, S., et al., 1999, ApJ, 517, 565 

Peterson, J.B., Bandura, K., Pen, U.-L., 2006, in proceed- 
ings of 41st Recontres de Moriond: Contents and Struc- 
ture of the Universe 

Platania, P., Burigana, C, Maino, D., et al. 2003, A&A, 
410, 847 

Prochaska J.X., Wolfe A.M., 2009, ApJ 696, 1543 

Reich, P., & Reich, W. 1988, A&ASS, 74, 7 

Riess, A. C, et al., 1998, ApJ, 116, 3 

Richards, E. A. 2000, ApJ, 533, 611 

Roberts, M. S., 1974, Science, 183, 371 

Rohlfs, K., & Wilson, T. L. 2004, Tools of radio astron- 
omy, 4th rev. and enl. ed., by K. Rohlfs and T.L. Wil- 
son. Berlin: Springer, 2004, 

Schlegel, D. et al., 2011, arXiv : 1106 . 1706 

Scott, D., & White, M. 1999, A&A, 346, 1 

Seiffert, M.; Mennella, A.; Burigana, C; Mandolesi, N.; 
Bersanelli, M.; Meinhold, P.; Lubin, P. 2002 A&A, 391, 
1185 



20 R.A. Battye et al 



Seiffert, M., Fixsen, D. J., Kogut, A., et al. 2011, ApJ, 734, 
6 

Seo, H.-J., & Eisenstein, D. J. 2003, ApJ, 598, 720 

Seo H.-J., Dodelson, S., Marriner. J., McGinnis, D., Steb- 

bins, A., Stoughton, C, Vallinotto, A. 2010, ApJ, 721, 

164 

Seymour, N., Dwelly, T., Moss, D., et al. 2008, MNRAS, 
386, 1695 

Springel, V. et al. 2005, Nature, 435, 629 
Szapudi, I., Prunet, S., & Colombi, S. 2001, ApJL, 561, 
Lll 

Tegmark, M. 1996, MNRAS, 280, 299 

Vernstrom, T., Scott, D., & Wall, J. V. 2011, MNRAS, 415, 
3641 

Weinberg D.H. et al. 2012, arXiv: 1201 . 2434 

White, R. L., Becker, R. H., Helfand, D. J., & Gregg, M. D. 

1997, ApJ, 475, 479 
Wilkinson, P. N., 1991, in Proceedings of the 131 st IAU 

Colloquium, ASP Conference Series, 19, 428 
Wolleben, M., Landecker, T. L., Reich, W., & Wielebinski, 

R. 2006, A&A, 448, 411 
Zwaan, M. A., Meyer, M. J., Staveley-Smith, L., Webster, 

R. L., 2005, MNRAS, 359, L30 



